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ABSTRACT 


A development of the in-plane open-loop rotational equations of 
motion for the proposed SCOLE in-orbit configuration is presented based 
on an Eulerian formulation. The mast is considered to be a flexible 
beam connected to the (rigid) Shuttle and the reflector. Frequencies 
and mode shapes are obtained for the mast vibrational appendage modes 
(assumed to be decoupled) for different boundary conditions based on 
continuum approaches and also preliminary results are obtained using 
a finite element representation of the mast-reflector sytem. The 
linearized rotational in-plane equation is characterized by periodic 
coefficients and open-loop system stability can be examined with an 
application of the Floquet theorem. Numerical results are presented 
to illustrate the potential instability associated with actuator time 
delays even for delays which represent only a small fraction of the 
natural period of oscillation of the modes contained in the open-loop 
model of the system. When plant and measurement noise effects are 
added to the previously designed deterministic model of the Hoop/Column 
system, it is seen that both the system transient and steady state 
performance are degraded. Mission requirements can be satisfied by 
appropriate assignment of cost function weighting elements and changes 
in the ratio of plant noise to measurement noise. 
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I. INTRODUCTION 


The present grant extends the research effort Initiated in previous 
grant years (May 1977 - Feb. 1984) and reported in Refs. 1-10 . Techni- 
ques for controlling both the attitude and shape of very large inherently 
flexible proposed future spacecraft -systems are being studied. Suggested 
applications of such large systems in orbit include: large scale multi- 

beam antenna systems; earth observation and resource sensing systems; 
orbitally based electronic mail transmission; and as in-orbit test models 
designed to compare the performance of flexible systems with that predicted 
based on computer simulations and/or on scale model Earth-based laboratory 
experiments. 

This report is subdivided into five chapters. Chapter II begins with 
a preliminary development of a two dimensional model of the rotational 
equations of motion for the proposed Spacecraft Control Laboratory Experi- 
ment - SC0LE^\ correcting some inconsistencies contained in a similar 
development appearing in last year's final grant report^. This development 
is based on the expansion of the Eulerian moment equations assuming that 
the Shuttle and the reflector are rigid bodies, and modelling the mast 
as a connecting flexible beam. Calculations are then performed to obtain 
the frequencies of the fundamental and subsequent bending and torsional 
modes as well as examples of the corresponding modal shape functions. 

Bending modal frequencies- are calculated based on different assumed boundary 
conditions: (1) where the flexihle mast is modelled as a cantilever 

attached to the Shuttle end such that the displacement and slope at the 
Shuttle attachment point are zero; (2) where the mast is considered as a 
flexible beam accounting for the rotational inertia of both the Shuttle 

^References cited in this report are listed separately at the end of each 
chapter. 
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and the reflector as end bodies. The cases of in-plane and out-of- 

» 

plane bending are treated as separate decoupled motions. Frequencies 

are also approximated through a preliminary application of the STRUDL 

finite element algorithm where the reflector as well as the mast are 

* 

assumed to be constructed of the same material. A comparison can then 
be made with the previous NASA results 11 and preliminary results pre- 
sented by the Harris Corp,. in the SCOLE Workshop 12 held at Langley 
during December 1984. Finally, Chapter II commences with a preliminary 
stability analysis of the open— loop ln-<plane SCOLE dynamics based on an 
application of Floquet's theorem. Analytical results can be obtained 
for two special cases: (1) where' the offset of the mast attachment 

point on the reflector is set to zero; and (2) where the gravity-gradient 
torques are not included, but the reflector attachment offset can be 
non-zero. 

In the following chapter the preliminary review of stability techni- 
» 

ques that can be applied when time delays are present in the implementation 
of control inputs, presented in last year's report, is now extended to 
include representative numerical results. . Examples considered include 
a second order controlled harmonic oscillator system and a fifth order 
system based on the dynamic model of the F-100 turbofan engine. These 
examples illustrate the potential instability which could result even 
for delays which represent only a small fraction of the period of natural 
oscillation of the various modes in the uncontrolled system. 

Chapter IV is based on a paper to be. presented at the Fifth VPI & 
SU/AIAA Symposium on Dynamics and Control of Large Structures and extends 
work previously initiated during the 1982-1983 grant year and partially 
supported during the Summer of 1983 on this grant. The evaluation of the 
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expected performance of the Hoop/Column antenna system subjected to 

stochastic Inputs is now extended to Include simulation of the steady 

state RMS errors in addition to the transient dynamics previously 

reported (Chapter VI of Ref. 10). The Kalman filter algorithm of the 
13 

ORACLS package is used to develop control laws and simulate the estimate 
of the state in an optimal LQG fashion. The results of Ref. 10 are also 
extended here to include the effects of non co-location of actuators and/or 
sensors. As a specific example the actuator (or the sensor) assumed to 
be mounted on the hoop assembly is then removed to examine the projected 
effects on the transient , steady state (RMS) errors , and the estimator 
performance. 

Chapter V describes the main general conclusions together with future 
recommendations. The effort described in Chapters II and III is being 

14 

continued during the 1985-86 grant period in accordance with our proposal 
and subsequent discussions. 
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II . MODELLING TECHNIQUES FOR THE SPACECRAFT CONTROL LABORATORY EXPERIMENT 


The transfer of large, massive payloads into Earth orbit is cur** 
rently accomplished with considerable propulsive and control effort. 

As a result, spacecraft designers must strive to minimize a large struc- 
ture's mass. Consequently, many of the future spacecraft will be very 
flexible and will require that their shape and orientation be controlled,^ 
The problem of controlling large, flexible space systems has been the 
subject of considerable research. Many approaches to control system 

synthesis have been evaluated using computer simulations including u pre- 

2 3 

limlnary synthesis of control laws for the proposed Hoop/Column System. * 

Ground experiments have also been used to validate system performance 

* 

under more realistic conditions but based on simple structures such 
as beams and plates.^ In a recent paper, SCOLE (Spacecraft Control Labora- 
tory Experiment), Lawrence W. Taylor Jr. and A.V. Balakrishnan described 
a proposed laboratory experiment based on a model of the Shuttle connected 

to a flexible beam with a reflecting grillage mounted at the end of the 
5 

beam (Fig. 2.1). The authors stressed the need to directly compare 
competing control design techniques, and discussed the feasibility of 
such direct comparison. Concern would be given to modelling order re- 
duction, fault management, stability, and dynamic system performance. 

With this paper^ as a background, the purpose of the study proposed 
here is to model the system in different phases where each successive 
phase would represent a mathematical model successively closer to that 
of the actual laboratory system. 
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It is anticipated that this (multi-year) study would consist of 
five parts, the first of which would consist of a literature survey 
during which the investigators would familiarize themselves with dif- 
ferent mathematical modelling techniques. 

During the second part, the system would be successively modelled 
as follows: 

a) The Space Shuttle as a rigid body; the reflector mast as a flexible 
beam type appendage; and the reflector as a rigid plate. The mast shape 
functions are actually solved from the fourth order non-linear flexural 
beam equation with different boundary conditions imposed on both the 
Shuttle and grillage ends, b) Here the Space Shuttle would be treated 
as a rigid body body; the composite appendage consisting of the flexible 
reflector mast and also the continuous rigid reflector (grillage) could 
be modelled using finite element techniques. Then the composite system 
dynamics can be modelled using the hybrid coordinate technique which 
involves sets of matrix equations describing the motion of the main vehi- 
cle as well as that of any attached appendages. It is anticipated that 
within the second part of this study these different mathematical models 
would be developed in a form suitable for numerical simulation. 

During the third part, each of these models could be directly com- 
pared with the model proposed in the SCOLE paper , beginning with the 
simulation of the open-loop system dynamics. The fourth part of the 
effort would consist of the control law synthesis when the model can 
be described by linear system dynamics - l.e>, in response to small per- 
turbations induced on the system about the nominal laboratory configura- 
tion and orientation, or after a major slewing maneuver, to remove the 
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remaining transients which exist in a neighborhood of the new equilibrium 

orientation. Such a construction of control laws will probably be based 

7 

on the ORACLS software package. Strategies would be developed to con- 
trol the shape and orientation of the beam/grillage. 

First the controllability of the system could be examined based 
on the graph theoretic techniques already employed for a similar analysis 
of the Hoop/Column system , for different combinations of numbers and 
locations of the actuators. Next, control laws can be constructed based 
on the techniques of optinial control theory, and studies can be performed 
comparing transient and control effort characteristics for a variety 
of system parameters and weighting matrix elements. 

Finally, the fifth part would focus on the slewing maneuvers to 
accurately point the reflector at a specific target *in a minimum lapse 
of time. For simple maneuvers (single axis) attempts would be made 
to analytically determine the slewing control law; for more general maneu- 
vers, numerical techniques would be implemented. 
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XI. A Development of the Two Dimensional Model - (Eulerian Moment 
Equations) 

The SCOLE system is assumed to be comprised of three main parts 
(Fig. 2.1): 

i) the Space Shuttle Orbiter with its center of mass located at 
point 0 ^; 

ii) the mast, treated as a 130 ft long beam, connected to the 
Shuttle at 0 2 and to the reflector at 0 3 f 

ill) the reflector, considered to be a flat plate with its center 
of mass at O 4 . 

The preliminary analysis presented here started before it was speci- 
fied 8 that the Interface point between the mast and the Shuttle is at 
Oj. Therefore, in what follows, a position vector appears which de- 
fines , where 0 2 is the assumed interface point. 

In the following analysis, the angular momentum, of the entire system 
is evaluated at point 0 ^ and the dynamics Include the lateral displacements 
of the beam. 

II. A. 1 Angular Momentum of the Shuttle with Respect to Point 0^ 
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Consider a point, P, of mass., dm, at an arbitrary position in the 
Shuttle such that O^P ** r. The elemental angular momentum of the mass, 
dm, is given by: 

c<Ho, s ^Opj^olm = r** ji.(R+rO| R dm 
a ^ R-fc + Rt^c 1/ ^ ( o3c- b) 3 x'r ] din (4.1) 

The total angular momentum for the Shuttle is obtained by integrating 

' •» 

Eq. (2.1) over the entire mass of Shuttle as: 

Ho, B il ' K) 

The first and second integrals appearing its the right side of Eq. (2.2) 


vanish because the center of mass of the Shuttle is at point 0^ , 


-h -V 

Since r • j =0, Eq. (2.2) takes the form: 


H V«. = («*" W $„r dm = Is 

where II S is the Inertia tensor of the Shuttle at point 0^ and = 


(w c - 0 ) j* 
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Consider here an element of the mast located at point, P^, with 
mass, dm. The elemental angular momentum of such an element is given 


hy; =jo,ft x £ q9,J r (dm 


if one notes that 


* ~ — S 

OP, r K + To+Cj 


(2.4) 


(2.5) 


then, Eq. (2,4) may be expanded according to: 

Hn/Oj z x J +jh)dx) ( 2 . 6 ) 


is expressed using the relationship between 


dffi 
dti&o 

the rate of change of a vector, w, in an inertial (R q ) and rotating 


(R.,) frames, i.e. 


(2.7) 


After substitution of Eq. (2.7) into Eq. (2.6) and integration term by 
term, one can develop: 

Hm/o, r MmftjJc A + 

L Cdn(uob t$) dcrsfit + C +J) cvd/i£ +8 

-f * tO /»« ( rtirt < p ) j 1(A +$ o < d \ ftt ) 

+ ! ( Ac»/}C + B /» tfl-C -A + cJj + lfa - P-vJjj 

AO /XYiUAdLcL - /{ CmjQxf d C CcA/fyC + 
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II. A. 3 Angular Momentum of the Rigid Reflector with Respect to 



0 | 

Let 0^ be the center of mass of the reflector., and 0^ the inter- 
face point between the reflector and the mast. The distance, X, between 
0^ and 0^ is constant since the reflector is assumed to he rigid, at least 
for this analysis. 

Let us now consider an element of mass, dm, of the reflector lo- 
cated at an arbitrary point, elemental angular momentum of that 

element of mass can be expressed as : 



•4* •'V 

0^2 and °2' P 2 can expressed as: 




( 2 . 11 ) 


Eq. (2.9) may be expanded according to 

1 


dH r / C/ =. [R,tOc* x) j| (R.+ Rii (x-tx.) t s J 


Once more, 
expressed using Eq. (2.7): 


d"h*3l"k+ S2 *y*** 


cttr 

After substitution of Eq. (2.7) into Eq. (2.11) and integration terra by 


term over the entire mass of the reflector, one arrives at 

7?r/ 0/ = / Mr RR t tO c (<X+ 9) t <A. v 

ft*, (ct+9) + (<*-»-«) £&r + e >)Jf 


( 2 . 12 ) 

A 

where is the moment of inertia of the reflector about the j axis 
taken at point 0^. 
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XI. B.l Moment Equation 


The angular momentum of the entire system about 0^ is obtained 
by summing the angular momentum of each part about Op i*e. 

4 



(2.13) 


(2.14) 


where N Is the sum of all the external torques, acting on the entire 
system, about an axis through point Op 

At this stage of the analysis, it is assumed that the center of 
mass of the Shuttle moves in a circular orbit, l.e. 


j n- 


(2.15) 


Taking into consideration the coincidence between points 0^ and 02» 
Eq. (2.14) is expanded using once more Eq. (2.7) and the following 


result is obtained: 




- §( I»i + 

_ 4 Mr(X l +R.t)) * («+*)k 4 Mm t*. 6 
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^ A jtM-ftfc - Sfi03 -t- C /Qwhjtyt 4 J> cedtjfyC 

_ Mr (&Ht <*.(*+*)+ XK<a)J 
+ Mm C& %!(<*+&) ^ fi*(cOb+(p J 

_ B cnftt ■+ C './*J> /%C +£ &U/& ■+ 8- d)J 

_j. M** CO yCcn Cioir + (f) foHj%C + £yfa*h/%£ - B coy$C 

4 2 ) ccvfijfiC ^ ^ / /) mfii •/ B _ £ an^feC 

& a 1 

™ — A •+ c J 


RftCb 


ex's 


$ =r 


X = 


JL On (boir+fi) 

-a? />*, (tat +<p)Wi) 
(u>t+(j>)<{(t) 


_ ccr 


MMn 


( 2 . 16 ) 
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II. B,2 Expression for q 

In the moment equation, Eq. (2.16), one notices integrals involving 

*+* 

q, the transverse displacement vector, and its first and and second deriva- 
tives with respect to time. It is, therefore, necessary to develop an ex- 
pression for q. 

II, B. 2.1 Relation between q(x,t) and y(x,t) 



Consider the beam in its deflected configuration*y(i,t) is the 
deflection of the reflector-end of the mast at an arbitrary tlme^ t; 
y(x,t), the deflection of an arbitrary point on the mast at the same 
time. 

A A 

From Fig. (2.1), k^ * k 2 = cos a (2.17) 

Assuming ct small, tan a can be expressed as 


;U« a SIS * Of = 

nr 


(2.18) 
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From Eq. (2.18) one derives 

CjCxjfc) = X _ y (z, b) 

•V 

or 

<yCx,t) = <vac- ^(x,t; 

Evaluation of y(x,t) 

Assuming separability of the variables, the beam equation, 

El. 91 _ n 

JA" 'J* 4 "M 1 

is solved to yield solutions of the form: 

£(x,fc) s |W $U.) 


f(t) « E sincut + F coswt with to ■ frequency of the vibration 
and <ji(x) » A cosgx + B singx + C coshfix + D sinhfJx 
When the following boundary conditions are assumed: 
a) y(0,t) » 0 ; b) y’(0,t) » 0 
c) El y MI (4,t) - -Mr y(4,t); d) El y M (&,t) = 0 


where 


i* - '' T.fr and u s 


Ur 


these can be expressed in the form: 


* A+ <5 8 so) 
y a + r& s oj 


M = 0 

T 6*0 
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where M f 

a * sin0Jl - sinh0& - 0 (coa0£ - cosh02.) 

M- 

& “ -cosflJl - cosh0£ - 0 (sin0£ - slnh02.) (2.26) 

Y “ cos0£ + cosh0& 
cr » sin0Jl + sinhBJl 



8 

For the SCOLE system, the following parameters have been supplied : 
pA / ■ 0.09556 slugs/ft 
El - 4. 0x10 7 lb-ft 2 
Mj. - (400/32.2) slugs 
% = 130 ft. 

For non-trlvial solutions for A and B, det C must vanish. The 
values of for which det C 3 0 are computed and substituted back into 
Eq. (2.28) to obtain the frequencies of the different vibrational modes 
(Table 2.1). 

The same values of 0 are substituted into <|>(x), (Eq. 2.23)^whlch is 
normalized with respect to its maximum value and the normalized mode 
shapes plotted (see Table 2.1 and Figs. 2.2 - 2.6). Note that the ranges 
of frequencies obtained in Table 2.1 are higher then those prevously pre- 
sented in the April 13, 1984 oral presentation due to previous inconsis- 
tencies in dimensional analysis of some physical units. 
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TABLE 2.1 


Values of 0 and Natural Frequencies (HZ) 
for the First 8 In-Plane. (Pitch) Bending Modes 


i 

1.874599 

4.6929 

7.8519 

10.997 

14.1309 

17.276 

20.4229 

23.555 


to (Hz) 
.677828 
4.245 
11.884 
23.3128 
38.4933 
57.5283 
80.4045 
106.958 




















XI. C Frequencies of the Lateral Vibrational Modes when the SCOLE 

System is Modelled as a Free-Free Beam with End Bodies having 
Inertia 

The solution to the beam equation (2.21) is again considered 
« 

and the following boundary conditions assumed: 

1. The shear force at either end is equal to the mass located at 
that end multiplied by the acceleration of the interface point at that 
end . 


This boundary condition combined with the equilibrium conditions 
yields 

- EX 9 s s Ml 

9z 3 ~ 

at the Shuttle end, 

EX . Ms s-Mtal'ufitjt)/ 

' ZZO 'Jt 2 ' fx *- 0 ^ f'Z-O 

at the reflector end 

El - Hr nf n - Mr to)* utZjt)/ (j.30j 

(xziio * 7 fr l kx/SQ ^ 

where CD 2 -- A 4 El 

1 fA> 
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2. Next, expressing the equality between the moment at an end point 
and the product of the Inertia of the mass at that end by the angular 
acceleration of the Interface point results in: 

J; OCzjtr) = El & 

CofltAJL QCXj t) - ) 

<~?Z‘ 

At the Shuttle end: x “ 0, this Is expressed as: 



the same boundary condition at the reflector end translates as: 

U> 2 lr QXfat) / - - ei ItMld h' (J.&&) 

Ox /zzho /xz-nc 


After performing the requir'd differentiation of the assumed 
solution of the beam equation (Ec . 2.23), one arrives at the following 
system of four equations with 4 unknowns, A,B,C, and D, 


Eq. (2.2 A + B + M&& C -2>sO C2.33 ) 

J’A 1 SA> 

Eq. (2.30) 
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Eq. (2.31) =^> - A + Is ft * 6 + & + D 2 0 (J. i$) 

J A' Sh } 


and Eq. (2.32) ^ 



~ ff. & l /*tod\ fit] C. -h ' ~ o (£• 


Equations (2.33) - (2.36) can be recast in the matrix format as 



For non-trivial solution of <(> (x) (Eq. 2.23) the determinant* of M(f3) 
must be zero. A computer program was written, and the values of (3, 
solutions ol : the nonlinear equation det [M(f3)] = 0, obtained. 

These values of 6 were substituted into 



to derive the frequencies of the inplane and out-of-plane lateral 
vibrational modes. The results are given in Tables (2.2) and (2.3). 
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Table 2.2 Values of 0 and Natural 


Frequencies (H ) for the first 9 
z 

In-plane (Pitch) Bending Modes 


3 


0.0097 

0.3065 

0.0310 

3.1308 

0.0549 

9.81922 

0.0789 

20.2809 

0.1030 

34.562 

0.1271 

52.6283 

0.1512 

74.4794 

0.1754 

100.229 

0.1995 

129.664 
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Table 2,3 Values of 3 and Natural 
Frequencies (H ) for the First 9 
Out-Plane (Roll) Bending Modes 


8 

“<V 

0.0103 

0.3456 

0.0310 

3.1308 

0.0549 

9.81922 

0.0789 

20.2809 

0.1030 

34.562 

0.1271 

52,6288 

0.1512 

74.4794 

0.1754 

100.229 

0.1995 

129.664 
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II. D. Derivation of the Frequencies of the Torsional Vibration, 
SCOLE Configuration. 


Assuming the mast to be a circular shaft, the torque at any point 
on the shaft is given by 

7 “_ qx ^ y&jb) 

^ x 

where G is the modulus of rigidity and I the polar moment of Inertia 
of the cross sectional area of the beam. This torque is opposed by the 
inertial torque y^ ( ty(X jfc)) 

* Qt 2 

where p is the density of the beam. For equilibrium, 

Gi qjbjt) _ j y -o 


Assuming the separability of the variables, Equation (2.38) is solved to 
yield, solutions of the form 


$(*,*) = ftt) (ptz.) 

f(t) = a cos (ait) + 8 sin (rnt) 

$(x) = A sin x + B cos ay'f T/^x 
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Boundary Conditions 


Writing that the torque, T, at either and of the beam equals the 
moment of Inertia times the angular acceleration of the interface 
point yields: 

GT - 2i Vj&t) O +o) 

Oz 

Equation (2.40) along with the equilibrium of the shaft gives: 
for the Shuttle end: x “ 0 

Ql j =: - Is LO 1 - u(x.,t)l (A.+O 

Qx* /XzO /Xzo 

for the reflector end: x = % " 130 

- GI - - 2r (x? z Q (Xj t) / (1*4* ) 

Qz /Zs/30 /XS./2C 

After substitution of equation (2.39) into equations (2.41) and 
(2.42), one arrives at: 

Eq. (2.41) a a i + b r s to*, o (a. 4-a) 


Eq. (2.42) *=> 
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Equations (2.43) and (2.44) can be recast in matrix format 



For non-trivlal solution of equation (2.39) one must insure that the 
determinant [P(w)] is equal to zero. 

The values of w for which det [P(oj) 3 a 0 correspond to the 
frequencies of the torsional vibration. A computer program was written 
to solve this deterrainontal equation and the frequencies for the tor- 
sional modes are listed in Table 2.4. 
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Table 2.4 Values of Natural 
Frequencies (H z ) for the First 9 
Torsional Vibration Modes 

0.0305 

39.99 

79.98 

119.97 

157.97 
199.96 
239.55 
279.94 
319.939 
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XX, E Preliminary Calculation of the SCOLE Appendage Frequencies 
based on Finite Element Techniques 

For this application both the reflector and the mast are assumed 
to be a single flexible body. This body is considered to be comprised 
of two types of elements; (1) beam elements; and (2) triangulur plate 
elements. The actual finite element model (FEM) is described ns follows; 
Mass distribution 

f 

Space Shuttle 6,366.46 slug 

Mast 12.42 slug 

Reflector 12.42 slug 

The masses of the reflector and the mast being so small (.39%) as compared 
with the mass of the orblter, which in this analysis is assumed rigid, 
the system could be modelled as a cantilever beam (mast) with a mass with 
inertia (reflector) at its end. Also, the reflector in this section is 
going to be assumed flat with a constant thickness small as compared with 
its characteristic dimensions. 

The dynamics analysis of the STRUDL software package, which uses a 
physical analysis to solve the equations of dynamic equilibrium, is used 
to generate the eigenvalues, the frequencies, and the periods of the system. 
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System Geometry (Model) 



The beam (mast) will be divided Into 3 beam elements (each of 
43.33ft length) having a mass of 4,14 slug to be lumped at the 
ends of the elements. 
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Coordinates of nodes for the System 










Results - Conclusions 


The following results have been obtained (Table 2.5). They show 

that the system is less stiff in this model as compared with previously 

developed NASA and Howard University continuum models and also that re- 

9 

cently described by the Harris Corporation. 


Table 2.5 - Modal Frequencies (H z ) 

Obtained by Implementing a FEM 
of the Preliminary Model of SCOLE 
(Poisson's ratio “ 0.3 assumed) 


0.157 

0.275 

0.782 

1.083 

1.232 

1.386 

80,09 

107.24 

107*24 

265.99 

421.50 


2.31 



II. F Linearization of the Equation of Motion-Floquet Analysis 

d d 

Let T the dlmenslonaless time be equal to W c t; ^ ; 

si * £_ 

2 2 
dx^ dr 

Dividing each term of Eq. (2.16) by M m £ yields 

- ® { 1 %~ i + X/ i + ^ ^ ] 

+ [(f - ^ |) tJt - ** £ <£) j ( p M>C u> fc ■* d>) - e +ij>)) 

- ft* ft§>«) 

+ flt + ( j4’» “ A** Cu3t+<(>) t»(cot+^) [^wJI 2 l^c 

~ !dr <*>!** Cwt + ij)} If <l) tlr x & tOc = o (2.+6J 

Mtn i l J Mwt t -t 

1 ^ t f +yk J 


ItH 

Let now C 


M_ 


where n “ ^ ~ X 

m 


% - (| - 

where ifi (&) =» A + B + C 603$. ^ H* /Smli 

m A ★ B Hr C/HwJfl ^fC •+ £> Cw& ^ *t B — 1 ) 

C., . ( :r ‘ r /M w ci + t / u.A*) l KO + 'h«)//4d* 

Where «p 2 (a)- 4 y&li - B CKy$t H- C/?tvilq ^ ■+ 2> cjp&fbt 

^3 00 = A 6 (n B />wi - C ecr^ - 1 ) y%£ — A 4 C 


2.32 



- (<n<)) {| j t -fy. 5 

- ykAfi 


and Cg " 


C*) /*- ^ £ 


Eq. (IM6) can be written as 

- C| S 4 ^(£ure^k+4?j) - 03 c C s ©■ 4 Ld l £± CcoQjOtid) 

ox ' • 

- CJ'.uJCj Intuit +$)Cjotu>b ity)- 00 (A: Cf Jutf ttrJj^rO 

C 2-4?) 

Introducing the dimensionalesa time T «* ai c t and dividing Eq. (2.47) 


by ( 1 )^ one arrives at 
c 


- OO^C, W + (4 1 C t L (&w»Cujfc+<|>J) - ai l c e & & | 

+ CjO Z C^ - UkMCq /fa fat entity )-<jQU)cC£l^foki$) = Oi 

a^si 


Introducing now the new parameter fl - — yields 

c 


- G, &" + fi. d [& to (ilt+£t>)] _ C 5 dgt 


op (£ 2^40 - 

4 


— *^1.^4 /hw, (XLT *■ 3 ^ Ii<— * Jja ) &^ s O 
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Parametric Study of the System 


Let us assume that the interface point between the reflector 
and the mast is at the center of mass of the reflector 
-> X ° 0 A = 0 = C 5 “Cg 

Under this assumption, the equation becomes 

Vs* £[* CO 027*4>)J + C* CEP 

C\ at c.| 

- fii - i(Ii-TiO$ J = O (^.S-o) 

AC i C>| 

which yields the following first integral 

^ 9 * + Ci +4>)J t j2 # £j /&w(.0x*<b) 

C| c r 

C/P ||z + (p) J| * ^ = ^ i.Z-’S’ 0 

C| Ci 

This equation can be plotted in the phase plane (0^6) for different 
values of p and Q. 


Floquet Analysis 


The angular motion about an axis perpendicular to the orbit plane 
is described by: 

&“= o.z j ® 

iZSZ) 
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This equation can be recast into the following matrix format 



( 2 -S3) 


where 


DPCOJ : 


-9$ + 9j- cjnD.Z 

C| C\ 


. [ & n fatfic) 4 i (r„ -r„)J 


i 


0 


Since P(t) is a matrix with periodic coefficients, the stability of the 
motion will be analyzed in what follows using the Floquet theorem. 


Ca3e__l. No gravity gradient, No offset 


-f>Cc) 


£*> cm£lX. 

Ci 


1 


[z(t)3 » CP(t)] [z(t)3 


- -Q.I>va O.Z 
c-i 

o 
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z., 


12 


-21 


■'22 


P 11 Z 11 * P 12 Z 21 
P 11 Z 12 * P 12 7 '22 


P 21 Z 11 + P 22 Z 21 


P 21 Z 12 * 


22 22 


( 1 ) 

( 2 ) 

(3) which becomes Z 2 j_ ■ Z^ since p 2 ^ 

1 and p 22 » 0 

(4) which becomes Z 22 « Z^ 2 


from (3) Z^ « Z^ substituted into 1 yields 

Z 21 " P 11 Z 21 + P 12 Z 21 

•• . 

Similarily from (4) Z 22 » Z^ 2 which substituted into (2) yields 
«• 

Z 22 " P 11 Z 22 + P 12 Z 22 
If one notices that P^2 ** Zt P 11 

* 

Chen Z 21 - p u Z 21 + P U Z 21 - ^ <P U Z 21 ) 


and Z 22 - P 11 z 2 2 + Pu Z 22 " dt <P 11 Z 22 J 


These two last equations are integrated and the following result 
for Z 21 z 22 obtained 

Z 21 “ P 11 Z 21 + K 1 

Z 22 * P 11 Z 22 + K 2 

4 

but from (3)^ Z 21 « Z.j^(t) and 

# 

from (4) ; Z 22 - z 12 (t) 

therefore, Z 21 <0) - Z 1;L (0) - 1 » P n (0) Z 21 (0) + K ± 
or for ■ 0 1 ■ since Z 2 ^(0) » 0 

+ k l » 1 

Z # 2 2<°) - Z 12 (0) = 0 =.p n (0) Z 22 (0) + K 2 


r\ o tr 


or for <J> “ 0 


— « - K 2 since Z^CO) ** 1 


Z 21 " P 11 Z 21 + 1 


Z 22 " P U Z 22 ~ —■ 
C, 


Solution of the linear first order equation 


dZ 


22 
d T 


* P 11 Z 22 “ “ C 1 


( 1 ) 


< 12 * 

The presence of '22_ and P-^2^22 ^ ec l uat i on suggests a pro- 

dx 

duct of the type <^(t)Z^(x) 

bu t d[^ u ) = U *«* + 4> 


dr 




( 2 ) 


Multiplying (1) by (j)(x) yields 




(3) 


which can become 


if one can find (the integrating factor) such that 

44 =-4fu 

dr 1 


( 4 ) 


( 5 ) 
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.-=? JU <f>(X) = Jf n c lx = f_ Cj cmSlT itc 


■=. ~ ^ fa i-ilt + K en- 
C|£ 


4 >Ct) = exj> ^ Aw ax -+ kJ 


from 


4 -C<P «*0 « - Sj. 4 GO 

dx <1| 


one arrives at 


$Zii - j"- |* 4 Ct) dr 


or 


£zz = 


_L (_ £». d)(c) dz 
Ob) J c, T 


«?»(£)= e«cp[^-A*jat#kj |-£i 


According to Taylor's series development of a function 


*“ “f t-a ' ini ] ■ «"{'-£■ 1 -j 


which is integrated term by term to give 


*» ■ - tr “ffe " im -* J( r -tr c ‘- (it ) 1 s -*'t)g4 


since <^3 (O) r 1 - " & 


c?zi Z,~* S2r £Xf f SlL Om ? ("^ Ql + C - Si T^j. + 

Ci l^lc, J [ Cz. Ci 2 ^ iy i 
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Solution of Z 21 “ Pj ^Z^i+l where p 11 


ooa At 


dr 


(j)^ X>i\ ■+ <j) d , 

dr dr 


omA ^ 

dr 


( 1 ) « ( 2 ) ^ 

dr 


« 

=> Jw. cjj'Ct) =. - -£i. Aw iir + K."* or 
<j>’ Ct) = *x/|o /hit, m ■+ k J j 

fc=j^dt 


from 


°?ii = J_ =- Awfl'C- K J ] 


According to Taylor’s series 


<^(j) = e*j=^£i_^*atje K - e K , i- Sit -> t 



which is integrated term by term to yield 

* “f 1 -M c ‘- 1 1‘ * - - - 

t 

Zz\ &) = O o 

<=<miCO = Cl) 

^0= ex r te ^ at ][t r ■&/ -r -jf^-o 

4* I* Avvv XXX J 

1 LC| it 


It can easily be verified that (jq) — 

e *tfe r ~ tl S = 1 

and finally 

<??21 — 



f'-tV-D’l'*-] 

(Cos SZI ^ CnC2t - £2-^ eX/J» 

/ c i fitoSlZ) 
I^Cj J 
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With the use of a computer program, the eigenvalues of the CZ(t)] 
matrix are computed for t « a period and their modulus compared with 
1 to determine the values of the parameters for which the system is 
stable. The results of such parametric study are shown In the 
following stability diagram, Fig. 2.7, The large number of unstable 
points in the parametric space (ft, MU) are thought to be attributed 
to the absence of the gravity-gradient torque in the model. Future 
plans call for the extension of the Floquet analysis for the cases 
where a non-zero reflector attachment offset is considered 4 and also 
where both a non-zero offset and the effects of gravity-gradient are 
Included. 
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to be multiplied 
by 103 


80,73 


57.92 



38.91 


♦ 


♦ 


23.72 


12.36 

4.97 

1.65 

0.29 



Fig. 2.7 Floquet Stability Diagram - SCOLE Coi 
No Gravity Gradient. 


• Unstable points 

* Stable points 

O possibly stable points 

A - jr- 0.0 
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III. STABILITY OF LARGE SPACE STRUCTURES WITH DELAYED CONTROL INPUT 

The dynamics of a large space structure can be represented 
by; 

X(t) - AX(t) + BU(t) (3.1) 

where 

X ■ 2nxl state vector 

A ■ 2nx2n system matrix 

B “ 2nxm control influence matrix 

U ■ mxl control input vector 

» 

n “ number of modes retained in the dynamic model 
The matrix A for a freely vibrating large space structure has all imaginary 
eigenvalues and, thus, the uncontrolled system is marginally 3table. 

The system given in (1) can be stabilized by a state variable feedback 
control law of the form 

U(t) = KX(t) (3.2) 

and X(t) - (A + BK) X(t) (3.3) 

such that the closed loop system matrix , (A+BK), has the required eigenvalues. 

In practice, due to the high dimensionality of the sta:e vector and 
the use of digital computers for evaluating control signals., there may 
arise a delay in the control input which can be mathematically modelled as; 

U(t) - KX(t-T) <3. 4) 

The consequent need to verify that the system described by 

X(t) ■= AX(t) + BKX(t-T) (3.5) 

is stable if the system described by equation(3.3) is stable is one of the 

topics of Investigation carried out during the 1984-85 grant period. 
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The system described by equation (5) is stable 3,1 if and only if 
-M(A) > | |BK| | (3.6) 

where 

pCA) ■ W A * +A) 

I I BK] I . >^ ax ((DK)*(BK)) 

A > (B1Q denote the conjugate transpose A and (BK), respectively 
and ^ max (A) ■ maximum eigenvalue of o\acrix A. 

Applying the result(3.6) to undamped large space stricture a (Rafflx A 
has all Imaginary eigenvalues or some zero eigvfttf&lues tfor^aspoutfing to 
the rigid body modes plus the imaginary sJgettVi^tjSes correspond tog to thu 
flexible modes), we have 
M(A) ® 0 
and | | BK | | > 0 

and, thus, the closed loop system becomes unstable. 

In reference 3.2, the analysis of time lag systems for stability is 
carried out with the assumption that the original system without these 
delays is stable. 

The large space structures pose a special problem in that the original 
systen is marginally stable and the verification of the closed loop system 
stability with delayed control input has to be carried out through numerical 
simulation. 

Two numerical simulation studies are conducted to determine the amount 
of time lag the system can tolerate without becoming unstable. 

CASE 1 : 

The system under consideration is an Jto.tmonic oscillator representing 
an isolated vibration of a large space structure at a specific natural mode 
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and Is given by: 


» 

ar 

0 

2 

1 ~ 


•*r 

+ 

-o - 

“ X 2" 


L-W 

0 _ 


- x 2 - 


_1 


U 


(3.7) 


The above system is stabilized with a control law of the form 


U 


Co * - 2^2 



(3,8) 


L X 2.J 

The numerical simulation is carried out for a control with time lag given 
by 


U(t) 


[0, -2 SwJ 


^(t-T) 

X 2 (t-T)^ 


(3.9) 


with the following numerical values: 
co => 6.0 rad/sec 

£ - 0.5 (3.10) 

x 1 (0) ■= 0.5 

x 2 (0) = 0.0 


The time response of x 2 (t) i9 plotted for t=0, T=0.085, x=0.1 secs 
(DT=0.05 secs is the numerical integration step time) in Figs. 3.1, 3.2 
3.3, respectively. The system became unstable for T=0.1 which is approxi- 
mately one tenth of the natural frequency of oscillation of the uncontrolled 
system. 


CASE 2 : 

The second system considered for numerical simulation is the dynamic 

3.3 

model of the F-100 turbofan engine * . The uncontrolled or open loop 
dynamics are stable. 
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(3.11) 


The system dynamics are mathematically described as: 

X(t) - AX(t) + BU(t) 

The controller is of the form 

U(t) » KX(t) (3.12) 

The matrices A,B,K and the eigenvalues of the matrices A and (A+BK) are 
given in Table 3.1. 

The numerical simulation is carried out for the controller of the 

form, 

U(t) “ KX(t-T) (3.13) 

with the initial conditions 

x 1 (0) « x 2 (0) => x 3 (0) « x^O) “ x 5 (0) “0.1 
for T“0, t:= 0.001 and t= 0.0099, and given in Figs 3,4, 3.5, and 3. 6 ; re- 
spectively. The system became unstable for T=0,0099 which is a very 
small fraction of the period of the highest damped frequency of vibration 
of the uncontrolled system. 

Conclusion 

As the stability criteria for systems with delayed state variable 
feedback (as reported in the literature) are found unsuitable to apply 
to large space structures, which are marginally stable without control, 
the numerical simulation is carried out to determine the, tolerable time 
delay without the closed loop system becoming unstable. It is observed 
that even very small delays can cause the closed loop system to be unstable, 
thus demanding the necessity of very robust controllers. 
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Table 3.1 SYSTEM MATRICES AND EIGEN VALUES 
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IV. EVALUATION OF PERFORMANCE CHARACTERISTICS FOR A SPACE ANTENNA 
SYSTEM SUBJECTED TO STOCHASTIC DISTURBANCES 


Abstract 

The problem of controlling a stochastic linear system representing 
the Hoop/Column dynamics by minimization of a quadratic Gaussian per- 
formance index, appropriately weighted in both the state variables as 
well as the control inputs, is considered. An optimal control law for 
the .finite element model of the Hoop/Column structural system without 
damping is realized by combination of the Kalman filter and linear 
feedback. A parametric study shows that suitable combinations of plant 
and sensor noise characteristics, and state weighting matrices can be 
selected to meet the mission RMS pointing requirements. The effect 
of removing the hoop-mounted actuator is to cause an increase in the 
RMS errors along with the increased control effort. An increase in 
least damped modal time constant is also noted, when the hoop-mounted 
actuator is removed. The effect of removing the hoop-mounted sensor 
is to cause an increase in the RMS errors along with the degradation 
in the estimator performance. However, removing the hoop-mounted 
actuator causes a greater degradation in the system RMS performance 
than removing the hoop-mounted sensor. 

I . Introduction 

Orbiting large flexible space structures have been considered for 
use in future large scale communications and other fields. As the size 
of the spacecraft system increases and the ratio of weight to area de- 
creases, flexibility considerations become very important. This is 
in contrast to small space structures which are assumed to be rigid. 

One such large flexible space structure which has been proposed for 
future space missions is the Hoop/Column Antenna System which is de- 
picted in Fig. 1 in its deployed configuration. 

The Hoop/Column system'*', contains a deployable (telescoping) mast 
system connected to the hoop by support cables under tension. The hoop 
contains 48 rigid sections to be deployed by motor drive units. 
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The desired shape of the RF reflective mesh is produced by a secondary 
drawing surface using control cables. The reflective mesh is connected 
to the hoop by quartz or graphite stringers. At one end of the mast 
the electronic feed assemblies are positioned, whereas at the other end 
are the principal solar arrays connected to the main bus based control. 

In order to achieve the required RF performance a pointing accuracy of 
+(0,03-0.10) degrees RMS and a surface accuracy of 12mm RMS will be re- 
quired. A finite element model of the Hoop/Column structural system 
without damping, is taken as the basis for the controls analysis. In 
order to reduce the dimensionality of the system, for computational ease, 
a modal transformation is carried out, In this case the truncated sys- 
tem contains thirteen modes, comprising the six rigid translational 
and rotational modes and the first seven flexible modes. 

The controls analysis of the Hoop/Column antenna system requires 
specification of the types of actuators and their locations and orienta- 
tions in the structure. For this study point thrusters and torquers 
are assumed to generate the required control forces and torques. Since 
we are considering thirteen modes, it is convenient here to choose a 
maximum of thirteen actuators in this analysis. Controllability and 
observability considerations of the Hoop/Column system based on the 
proposed location of actuators as shown in Fig. 2 have been established 
using graph theoretic techniqes. 

An optimal control law is realized by combination of the Kalman 
filter and linear feed-back techniques (Fig. 3). The controls analysis 
is carried out assuming collocated and noncollocated sensors and actua- 
tors. The plant and sensor noises are assumed to be uncorrelated, zero - 
mean white Gaussian processes. 


II. Mathematical Formulation of the Problem 

The dynamic model of the Hoop/Column structural system in the ab- 
sence of damping can be represented as^ , 

M 

MZ + KZ - F„ (1) 

• c 

where 

M - 6n x 6n mass /inertia matrix 

K - 6n x 6n stiffness matrix 

Z - 6n x 1 matrix consisting of the displace- 
ments and rotations at the nodal points 

F - 6n x 1 control vector, 
c 

F ° B U (2) 

c c 
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where 

B„ - control influence matrix of order 6n x p 
c 

for 

p - number of actuators 

U - p x 1 matrix associated with the control vector. 

In the present model represented by Eq. (1) the number of nodes 
is equal to 112 (i.e. n * 112), corresponding to the number of nodal 
(grid) points in the FEM output and the number of actuators is equal 
to 13 (or 12) . 

To decrease the dimensionality of the system a modal transfor- 
mation is carried out defining 

Z « 4>q (3) 


where 4> is the matrix containing the eigenvectors of Eq. (1) and 
is of order (6n x m) , for m number of modes and q is a modal vector of 
order (m x 1) . In this case, m ■ 13. 

After using the transformation, Eq. (3), in Eq. (1) there results 

t 

* T M$ q + 4 |T K4’q » 4> T F c (4) 

The left hand side of Eq. (4) can be rewritten t , using the properties 
of the eigenvalues and associated eigenvectors as 


[ s m. ] q + L" K, ] q - 4> T F 
X — *■ c 


(5) 


where 




T 

$ X K$ 


diag C “ m ± ^ ], i » 1,:>..13 

diag [ 3 « [' ] f i - 1,2. .13 

Eq. (5) can be rewritten in the form 
X » AX + BU 

where 


( 6 ) 


X - 


q l‘ 

L q 2j 




:] 
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and 


i m ± r 1 * t b c 


The plant noise, w, may be incorporated Into Eq * ( 6) to model the 
stochastic linear dynamic system of the form: 

X • AX + BU + Gw (7) 

The measurements, Y , are assumed to be related to the state through 
the observation matrix, H, and the measurement noise, v, by 

Y - HX + v (8) 

Eqs. (7), (8), along with the cost function 

J «* lim {E CX T (t.)q.X(t f )3 + A f (X T QX + U T RU)dt]} (9) 

tf+c. f f f t 0 

completely define the stochastic problem,^ The minimization of the 
cost function yields the optimal control vector, U, obtained from 


4k 

U - -C X (10) 

A 

where X is the state estimate, and 

C » r" 1 b t k 1 (H) 

where is the steady state solution of the matrix Riccati differential 
equation, 

-Kj ■ K X A + a \ ± -K 1 BR“ 1 B T K 1 + Q (12) 

A C 

The state estimate X, is obtained, from 

X “ AX + BU + ]? (Y-HX) (13) 

with the filter gain, F, expressed as 

F “ PH T V" L (14) 

where, ? is the steady state solution of the filter matrix Riccati 
equation: 

P » AP + PA T - PhV’ 1 HP + GWG T (15) 


with 

V 5 (t-T) ■* E[w(t)w T (T)] 
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and 

V 6 (t-T) « E[v(t)v T (T)3 * (17) 

Substitution of Eq. (10) into Eq. (7) and Eqs. (8) atld (10) into 
Eq. (13) will yield the following differential equations which could g 
be used for simulation of stochastic optimal control systems (Fig. 3). 

X - AX -BCX + Gw (18) 

X - (A - FH - BC)X + FHX + Fv (19) 

IIA Simulation of the Steady State RMS State Components 

» 

After subtracting Eq, (18) from Eq. (19) a differential equation 
for the error results as follows: 

* e ■ (A - FH) e + Fv - Gw (20) 

Furthermore^ Eq, (19) can be rewritten as 

X * (A - BC) X - FHe + Fv (21) 

The covariance of the reconstruction error, P “ ECee^], and the co- 
variance of the state estimate, X,“ ELXX t 3, can be obtained from Eqs. (20) 
and (21) as, 


• T T -1 T 

P - AP + PA - PH v HP + GWG 


(15) 

• 

X- - (A - BC^ + X x (A - BC) T + 

T 

FVF 

(22) 

The covariance of X is given by 

ECXX T ] » Et(X - e) (X - e) T ] - 

*1 + P 

(23) 


Fgr the case when P ■+ 0 and X, + 0, the steady state variances of 
e and X and, hence, of X can be found. 

III. Possible Arrangement of Actuators for the Hoop/Column System 


Out of the maximum of thirteen actuators selected twelve actuators 
consisting of point actuators and a torquer are assumed to be located 
at selected grid positions in the feed assembly as well as on the column, 
while the remaining actuator is assumed to be. a point actuator mounted 
on one of the rigid links of the hoop assembly with thrust direction 
tangential to the hoop circle. 
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Table 1 indicates the relationship between the actuators and the modes 
they directly influence, the frequency of each mode, the generalised 
mass and the generalized stiffness. 


Table 1 - Relation between Actuators and Modes Directly Influenced 


Actuator No. 
(circled in 
Fig. 2) 

1,2,3, A 4 

Mode Affected 
( ) 

Feed Mast 
Torsion (12) 

Frequency 

Hz 

0.88976 

Generalized 

mass 

m i 

(lb-8ec"/in) 

723.522 

Generalized 

stiffness 

k i 

(Ib/in) 

22612.9 

5 

First Bending 
(about Y axis) 
(8) 

0.214246 

5.233 

9,h83 

6 

First Banding 
(about X axis) 

(9) 

0.270956 

3.073 

8 . 907 

7 

Surface Tor- • 
slon (10) 

0.506323 

0.305 

3.083 

8(Torquer) 

Yaw (rot. 
about Z axis) 
and First 
Torsion (7) 

0.0 and 
0.118835 

8.419 & 
153.157 

0.0 S 
85.385 

9 

Transl. along 
X axis A 
Second Mast 
Bending (11) 

0.0 and 
0.728873 

16.444 A 
1.993 

0.0 A 
40.887 

10 

Transl. along 
Y axis & 
Second 

Mast Bending 
(13) 

0.0 and 
0.919231 

8.925 A 
0.658 

0.0 A 
21.954 

11 

Transl. along 
Z axis 

0.0 

7.349 

0.0 

12 

Pitch (rot. 
about Y axis) 

0.0 

2.941 

0.0 

13 

Roll (rot. 
about X axis) 

0.0 

9.704 

0.0 
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NOTE: If there is more than one numerical value ifi any oolumn, the 

first one corresponds to a rigid mode while etfcpiut' corresponds to 
the flexible mode, respectively. 

IV. Numerical Simulation and Synthesis of CiiUtrol Law 

"“ " 1 1 ‘ " 1 ' it*- i — .1^^— 

For the proposed arrangement of actuators in the structure, a 
parametric study was performed showing the effect of varying the 
state penalty matrix Q (from 1001 to 100001) and fctofc A^ntrol penalty 
matrix, R(from I to 1001) on the least damped mtvs tJ? the system 
(Fig. 4). It has been concluded that Q “ 10001, itV a suitable 
design point from the stand point of minimizing the least damped modal 
time constant and maintaining a reasonable control effort.^ Further, 
the stability of the closed-loop system, consisting of the plant 
and the estimator, has been established for all combinations of 
parameters considered in this study, ? 

Plant noise covariances ranging from 0.00001 1(26 x 26) to 
0.0000001 1(26x26) and sensor noise covariances ranging from 
0.0000025 1(13 x13) to 0.00025 1(13 x 13) are assumed in the analy- 
sis. Selection of the range of plant noise covariances in the modal 
coordinates is based upon plant noise intensity ranging from 1 x 10”^ 
to 1 x 10“8 (lb-m). 2 Selection of the range of sensor noise covariances 
in the modal coordinates is based upon sensor noise intensities in 
the range of 1 x 10"5 to 1 x 10“? [(rad)2 or (rad/sec) etc].® Based 
on the accompanying analysis, it is seen that the estimated modal 
coordinates closely follow the actual modal coordinates, thus ensuring 
a good estimation process. ? * 

In this study, initial displacements of 0.01 are assumed in all 
the modal coordinates, .which correspond to the expected maximum allow- 
able perturbations in the linear range from the nominal operating re- 
quired RMS displacements; these amplitudes for the modal displace- 
ments are obtained through calculation of Eq. (3). 


Numerical studies show that when noise is included in the model 
that increasing the elements of the state weighting matrix by an order 
of magnitude increases the control effort required by an order of 
magnitude, but does not cause a significant improvement in the system 
response,' (But, it was found in the deterministic case that in- 
creasing the elements of the state weighting matrix causes a significant 
improvement in the system response. 3) 


A further parametric study was conducted to determine the range 
of system parameters which meet the mission RMS pointing accuracy 
requirements. Fig, 5 shows the relationship between the steady 
state RMS pointing errors (in some of the modal coordinates) and 
the state weighting matrices, for initially assumed noise covariances 
of w = 0.00001 1(26 x 26) and V = 0.0000025 I (13 x 13). 


4.7 



S. ANANTHAKRISHNaN, P.K. BAINUM, AND A.S.S.R. REDDY 


Fig. 5 shows a similar relationship for plant and sensor noise co- 
variances of W “ 0.0000001 1(26 x 26) and V - 0,000251(13 x 13), 
respectively. A comparison of Fig. 5 with Fig. 6 reveals that there 
is a decrease in the magnitude of the RMS errors in the latter case, 
for all values considered for the steady state weighting matrix ele- 
ments. In addition, the steady statn. RMS errors, for almost all nu- 
merical combination'' of the elements of the state weighting matrix 
and of the noise covariance matrix considered in Fig. 6, are within 
'ihe allowable maximum as specified by the mission requirements; as 
a contrast some of the RMS errors shown in Fig, 5 would actually 
exceed the mission specifications, [it is also important to note 
that the range of maximum actuator force amplitudes is about the same 
in both the cases,'] 

3 

In a recent related study , it was found that the hoop-mounted 
actuator play3 an important role in controlling the system. In this 
analysis, the effect of removing the hoop-mounted actuator (whose 
thrust direction is assumed to be tangential to -the hoop circle) is 
Kn cause a deterioration in the least damped modal time constants. 

The range of the maximum force amplitudes required, in response to 
initial displacements of 0.01 in all the modes, is also increased 
(Table 2). Further, amplitudes of the RMS pointing errors are increased 
substantially in all the coordinates as compared to the case where all 
12 actuators are present, as can be seen by comparing Fig. 7 with Fig. 5. 

The effect of removing the hoop-mounted sensor is to cause a 
deterioration in the least damped time constant associated with the 
estimator poles, thus degrading the estimator dynamics. There appears 
to be a redistribution and a smell increase :Ln the RMS errors as is 
evident by comparing Fig. 8 wit);} Fig. £. 

In another recent related study^, only torque actuators located on 
the mast are considered for controlling the antenna attitude and flexi- 
ble motions (using LQG theory as well as other techniques). However, 
it was suggested that hoop-mounted control devices might be effective 
in controlling certain torsional modes, which is made clear in our study. 
It should be kept in mind that the problems related to the design of con- 
trol hardware remain yet to be answered, 


V. Conclusions 


Control system synthesis is considered here for a large space an- 
tenna system and is based on stochastic linear optimal control techni- 
ques. From this preliminary analysis, it appears that the performance 
requirements can be met by using linear quadratic Gaussian techniques. 
Parametric studies show that suitable combinations of plant and sensor 
noise characteristics, and state weighting matrices can be selected to 
meet the mission RMS pointing requirement: s; here a definite trade-off 
exists between the increased complexity, cost, weight, and reliability 
of the system, and the possible gain in the system performance. 
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Table 2 



Comparison of Maxim im Actuator 
Q « 10301, R-I, q ± (0)-0.0X, 

Force Amplitudes 
i a l f 2 « • • 1 13 

W » 0 

13 Modes 
.00001, V » 

0.0000025 


Maximum actuator force 

13 Actuators 


’ “ Actuators 

amplitudes (pounds) 

13 Sensors 


13 Sensors 

f l 

0.3330 


2.3548 

f 2 

0.0570 


0.0453 

f 3 

0.2681 


2.2954 

f 4 

0.0570 


0.0453 

f 5 

1.3028 


1.3024 

f 6 

0.2865 


0.5520 

f 7 

1.2310 



f 8 (in-lb) 

0.0140 


0.0644 

f 9 

0.2859 


0.5328 

f 10 

0.1574 


0.1694 

f ll 

0.4086 


0.3693 

f 12 

0.3521 


0.3521 

f 13 

0.0660 


0.0993 
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The effect of removing the hoop-mounted actuator is to cause an. 
increase in the RMS errors along with the increased control effort. 

An in, crease in the least damped modal time constant is also noted, 
when the hoop-mounted actuator is removed, [it is interesting to 
note that removing the hoop-mounted actuator resulted in an increased 
control effort and also increased the least damped modal time constant, 
in the deterministic study also. 3 ] The effect of removing the hoop 
mounted sensor is to cause an increase in the RMS errors along with 
the degradation in the estimator performance. 

In order to arrive at more complete conclusions, factors such as 
time delay, and nonlinearities in the plant and sensors should be con- 
sidered. Since many of the current (and proposed) sensor systems pro- 
vide data in a discontlnous (discretized) format, another suggestion 
is to reformulate the current analysis in a discretized time basis. 
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First Bending 



Fig. 2. Proposed 'Arrangement of Actuators-Hoop/Column System 
(uncircled numbers identify grid points from finite 
element analysis) . 
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.agonal; Fig. 6. Steady State RMS Errors In Selected Modal 

Steady State RMS Errors In Selected Coordinates for Collocated Actuators and 

Modal Coordinates for Collocated Sensors - 13 Actuators/13 Sensors/ 13 

Actuators and Sensors 13 Actuators/ Modes - Effect of Further Decreasing W 

13 Sensors/ 13 Modes. and Increasing V . 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A two dimensional model of the SCOLE configuration open-loop 
dynamics has been developed and will provide the basis for the modelling 
of the three dimensional dynamics wich has already been initiated. 
Calculation of open-loop appendage frequencies is seen to be dependent 
on the boundary conditions assumed at both the Shuttle and the reflector 
ends of the mast, and also somewhat dependent on whether a (decoupled) 
continuum or (coupled) finite el ament model is used to model the append- 
age modes. The Floquet stability analysis, initiated during this report- 
ing period, should provide considerable Insight into possible parametri- 
cally induced open-loop instabilities which should be understood prior 
to the design of the control system. 

Numerical examples illustrate the potential instabilities that could 
result for large ordered systems (typical of LSST) even with input delays 
that are only a small fraction of the system's frundamental open-loop 
period. Such delays could be associated with actuators or other system 
hardware. This work should be extended to analyze the effects of input 
delays which could be associated with the SCOLE in-orbit or laboratory 
test scale model configurations. 

A parametric study of a stochastic closed-loop linear system dynamic 
model of the Hoop/Column system shows that suitable combinations of plant 
and sensor noise characteristics and LQG state and control weighting 
matrix elements can be selected to meet the mission RMS pointing and antenna 
shape accuracy requirements. 
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The important role of a hoop-mounted actuator in reducing RMS errors as 
well as improving the transient response characteristics, especially 
in the torsional type modes should be noted. Failure of the hoop-mounted 
sensor could result in an increase in RMS errors and also a degradation 
in the estimator performance. 
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